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We construct models describing interaction between a spin s and a single bosonic mode using a 
quantum inverse scattering procedure. The boundary conditions are generically twisted by generic 
matrices with both diagonal and off-diagonal entries. The exact solution is obtained by mapping 
O ' the transfer matrix of the spin-boson system to an auxiliary problem of a spin-j coupled to the 

spin-s with general twist of the boundary condition. The corresponding auxiliary transfer matrix 
■ is diagonalized by a variation of the method of Q-matrices of Baxter. The exact solution of our 

problem is obtained applying certain large- j limit to su(2)j, transforming it into the bosonic algebra. 

Introduction Models representing interactions between a single bosonic mode and spin degrees of freedom find 
i— T application in many different contexts. In atomic physics they describe atoms interacting with electromagnetic 
field [1] and many phenomena like spontaneous emissions in cavity [2] and Rabi oscillations [3] in two-level atoms 
are captured by an important representative of the models mentioned above, the Jaynes-Cummings (J-C) model [4]. 
. Recently the J-C dynamics was intensively studied in the research field of ions in harmonic traps [5] and then in 
quantum computation [6]. Finally this kind of models found applications in quasi-2D semiconductors in transverse 
magnetic field [7]. 

Toy model Hamiltonians representing a single bosonic mode interacting with a spin s are of the following type 

~0 H = H s ^ ph + 7 (S + a ! + S~a) (1) 

where H s - V h — uja^a + aS z + (3 (S + a + S a'); operators and a are bosonic operators that commute with spin 
i Q | operators S v , v — {z, ±}. The model (1) with 7 = (3 was originally proposed for quantum optics purposes to describe 
a dipole-like interaction in atoms-radiation systems and it is known as the Tavis-Cummings (T-C) model [8]; the 
£NJ , model reduces to the J-C one for s = 1/2. In solid state physics models of type (1) can describe certain quantum 
> ■ circuits [10,11]. 

There are two cases in which the model can be solved analytically: i)ln the limit s — > 00 there are exact results [12]. 
They were applied to study the entanglement across the quantum phase transition between normal to super-radiant 
phase [13]. ii) For "single photon" interactions the model can be simplified employing the Rotating Wave Approx- 
imation (RWA) that neglects the so called "counter rotating" terms: S + aJ, S~~a. Within the RWA the model (1) 
can be solved exactly [4,8,9]. For generic parameters a, ft, 7, to, and for finite s the model, as it stands in Eq. (1), is 
non-integrable. Merging the model into the main stream of the Quantum Inverse Scattering (QIS) method [14] con- 
stitues often a guide to discover unsuspected exactly solvable models with sufficiently generic interaction. According 
to this method the Hamiltonian is obtained as output of the procedure that remarkably ensures the integrability of 
the theory. In the simplest cases the QIS method provide integrable Hamiltonians after periodic boundary condition 
are imposed. The variety of the integrable models can be considerably enriched by considering more general boundary 
conditions [15]. By this is meant that the monodromy matrix is multiplied by non-trivial matrices that ultimately 
cause the presence of boundary terms in the Hamiltonian. Of interest in the present paper is the case of constant 
boundary matrix; this realize the, so called, twisted boundary conditions. For the type of models under consideration 
the QIS method was employed in Ref. [16-19] where nonlinear generalizations of H s - p h were studied. These gener- 
alizations were obtained by twisting the boundary conditions. The twist matrices were chosen as the same diagonal 
matrix for both the bosonic and spin degrees of freedom; finally in a certain sense (specified below) they are classical. 
As a result, although nonlinear, these models contain the standard interaction. Here a more general interaction is 
obtained applying more general boundary conditions: the twist matrices Kb and K$ are non-diagonal, different for 
the boson and the spin, and of "quantum" nature (see (5)) [20]. The Hamiltonian we found is Eq. (8). The QIS 
method pave the way towards the exact solution of the theory through Bcthc Ansatz. In the cases where there exist 
an obvious "reference" state a direct (algebraic) Bethe Ansatz approach can be applied. For the model we found here, 
however, there is no simple vacuum state since (8) does not commute with S z + a^a. A standard route to attack the 
exact solutions of such kind of spectral problems is to apply the technique that Baxter [21] invented to obtain the 
eigenvalues without knowledge of the eigenstates. We obtain the eigenvalues (the calculation of the exact eigenstates 
will be the object of a future publication) in the following way. We first define an auxiliary problem consisting of two 
spins with two distinct representations j and s; the boundary conditions are generically twisted; the spin j is affected 
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by an "impurity" v. We diagonalize the auxiliary problem by adapting the Baxter method to it. Then the solution 
of the spin-bosonic problem is obtained performing certain j — > oo limit (see (13)) in the results for the auxiliary 
spin-spin problem. The eigenvalues are given in (21) and the parameters Xk are fixed by (22). 

The paper is laid out as follows. In the next section we derive the integrable model. In the section III the exact 
solution is obtained. The section IV is devoted to our conclusions. 

Integr -ability. The starting point of the QIS method is to define quantum Lax matrices L(X) and a scattering 
matrix R(X) satisfying the Yang Baxter (YB) equation: R(X — n)L(X) £g> L(fij = L(fij ® L(X)R(X — jj,) , where A is 
the spectral parameter. For the present case, the Lax operators L we consider [17] are 



Lb(X) ■■'("-"-I 1 - 4-> ) • (3) 

each satisfying the YB equation with: R(X; if) = »7ll<g>ll+AP , where r\ £ IR and P is the permutation: PA<E)BP = B®A. 
The monodromy matrix is 

T(A) = K B L B {X)K S L S {X) , (4) 

where Kb and K$ are C-numbcr matrices that produce boundary terms (without "internal dynamics" , the matrices 
K not depending on A). Notice that we have two different boundaries each for the spin and for the boson. In 
Refs. [17,18] Kb = Ks is assumed; the T-C model (without counter rotating terms) H s - P h is obtained for r/ — ► 0. 
The matrix T(A) fulfills the YB relation: R(X - /-t)T(A) <g> T(n) = T(/z) <g> T(X)R(X - /x) , due to the fact that 
[R, Kb <8> Kb] = [R, Ks ® Ks] = holds for any numeric matrix because of the sl(2) symmetry of the ii-matrix. The 
transfer matrix is defined as t(X) :— trr Q \T(X) where trv ) means trace in the auxiliary space. t(X) is a generating 
functional of integrals of motion since: [t(X),t(fi)] = 0. For the present case the transfer matrix can be chosen as a 
polynomial in 77: t(X) — X^=- s rfhW'i tncn tnc coefficients of [t(X), t(/z)] = J2i=-2g V l ^l(\ A 4 ) vanish at any 77-power. 
The following assumption is crucial for our purposes: The entries of the matrices K depend on 77 

Kx a = K% + K^ ijV + ... , X = {B,S}. (5) 

The parameter 77 is usually called "quantum parameter" since it controls the limit how to recover the classical scattering 
matrix r(A) from the matrix R(X). In this sense our matrices K in (5) describe "quantum systems" which we couple 
to the the boson and to the spin at the boundary (twist matrices that are independent on r\ might be considered 
as classical boundaries). In brief, the main idea of our procedure is to play with the boundaries K in such a way 
that Ci(X,fi) = O [ti( A), ti(jj,)] = for certain I. In order to obtain the model we are interested in, the degree 
and the coefficients of the polynomials are fixed such that: i) t m (X) describes an integrable model, then tj(A) must 
be C-numbers for all I < m; ii) the model results containing the counter-rotating terms; Hi) the obtained operator 
is Hermitian. All these conditions translate in a system of equations for the entries of the matrices K; we shall see 
that these parameters will be the coupling constants of the Hamiltonian. It turns out to be sufficient to consider the 
entries of the K matrices to be linear in r\. Such entries are restricted to 

'(0) 7^(0) 



K (0) _ K (0) K {0) _ K (0) K (0) _ K (0) K {0) _ K (0) „(0) _ _^B21^_S22 / fi x 
"Bll — "B22 i "Sll — "522 ' JX B21 — "S12 ' "512 — ^ S 21 i- fv S12 ~ (0) ' W 



jfW _ *B22*S22 , A Bll /y(l) _ K W \ _ flSlL( K W K W ) I ^S22 _ K W 

^Sll— (0) + „(0) *--"-S21 ^SUJ (0) \^B12 tx B2l)^ ( ) ^ B 11 ^S2 

^B\\ B 12 "-B12 

We take t\(X) as the Hamiltonian 



h(X) = H = W{X)S Z + X(U + Vya + 2 Y + VUV(A - A) S x + ZX{a + a*) (7) 
-U(aS+ + tfS-) + V(a f S+ + aS~) - 2VlJV(a + a^)S z , (8) 

where the couplings are 
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W(X) = +(V- U)(A A - -j=) , (9) 

U = -Kf 22 Kf 22 , V = ~Kf 22 {Kf 2l f/Kf 22 (10) 

^ = ^B21 ("^B ll^S2l/^B22 — ^S2l) — ^B12-^S22 — ^B 22^5 12 ' (H) 

Z = Y + Kf 22 (Kf X2 + K&) + 2K^ 22 Kf 2l . (12) 

The coupling constants obey (9) for the model to be integrable; parameters A, A, X, Z can be set freely; the 
quantity UV must be positive. Nevertheless the rotating and counter-rotating terms can be adjusted to have the 
same sign by acting on the operators: S~ — ► S~e l % and a — ► ae l %; the third, fourth and last terms in Eq. (8) arc 
transformed accordingly. We observe that the simultaneous presence of rotating and counter-rotating terms preserves 
the integrability only if a further term (a + a^S 2 appears in the model (the term a + can be transformed out by 
a translation: a — > a + £ with £ = — Z/(U + V); the coefficients of S x and S z are shifted by £({/ — V) and — 2£\/UV 
respectively). Restricting the boundary conditions: Kb = K$ (non-diagonal) induces the further constraint U = V. 
In this case our model reduces to the Rashba Hamiltonian in a constant magnetic field [22] (the bosonic number 
labeling the Landau levels; see also Ref. [23]). 

The constants of motion of (8) are ti(0) and d\ti(X) (only two of H, ti(0), d\ti(\) are independent). d\ti(X) 
can be easily diagonalized : Ud\ti(X)U~ 1 \(f>) — ((U + V)(n + m) — MJ)\4>). In the \<j>) basis the Hilbert space of 
the Hamiltonian blocks into invariant subspaces labelled by the bosonic number n > and with S z \cf)) — m\<f>)\ 
M- = Z 2 /(U + V). This will be used to classify the excitations in the Bethe equations (22). 

Exact eigenvalues. To diagonalize the model (8) we define an auxiliary inhomogeneous spin problem. We use the 
property that a spin j-su(2) can be contracted to the Weyl-Heisenberg algebra [24] through the singular limit e — > oo 
of a Dyson-Maleev transformation 

-vJ~, ~2 J+ ^ JZ | ^ | flt ' a ' - fttft - yj ; ( 13 ) 

such a limit corresponds to j = —e 2 /2 — > oo. The bosonic Lax matrix is thus expressed as limit of a spin-j Lax 
matrix: 

L B (X) = lim ( — -)K(e)a y a z Lj(X - v)a y , (14) 

e— >oo tlS 

where K{e) — diag {r]s, {r/e)^ 1 } and 

the "inhomogeneity" parameter being set to 

v = -r)e 2 /2 + n^ 1 + A . (16) 

Thus the monodromy matrix Eq (4) can be written as T(A) = lim e ^ 00 T a (A) where T a is an auxiliary monodromy 
matrix defined as 

T a = KjLj(X - v)(jyK s L s (X) , (17) 

with Kj = —l/(7]e)KBK(e)a y a z . t a — tro{T a } can be diagonalized adapting the Baxter method [21] for off-diagonal 
twisted spin-s chain [25]. In the present case the "chain" consists of only two sites; the twist matrices are distinct 
and containing both diagonal and off-diagonal entries. The details of the calculations will be reported elsewhere. The 
Baxter equation reads 

t a (X)Q(X)=R™ ( A - sr? )(A-z/-jr ? )Q(A + r ? )+i?( ± )(A + S r ? )(A-z/ + ir ? )Q(A-r ? ) , (18) 



where = \ 



tr(Kja y K s ) ± ^{tr{Kja y K s )f +Adct{KjK s ) 

matrices and constructed in the standard way [21,25]; they fulfill 
eigenvalue r a {X)^ of the transfer matrix 



. The quantities Q(X) are (2s+2j+2) x (2s+2j+2) 
Q(A), Q(ji)] = {Q(X),t a (X)] = 0. Eq (18) fixes the 
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r a (A)W = R™ (A - s V )(X v j V ) 2 f{ X X ' +V 



A-A,- 

(A + s V )(X v + jr?) J] X ~^~ V , (19) 
i=i A Ai 

where the variables Xj are solutions of the equations 

(Afc - gr?)(Afc - t/ - jg) _ ^ A fc - A, - r? 
J?(±)(A fe + S r ? )(A fe -^ + j7 ? )- 11 X k -X t + ^ k-1,... ,2(8 + 3). W 

To obtain the solution of the bosonic problem we perform the algebraic contraction at the level of Eqs. (20). This 
can be done expanding R^ = Y^ m ^2m /( £7 ?) 2m an d taking into account of Eq. (16); a generalizations of the Bethe 
equations found in [17] arc obtained. By a further expansion i?2m = ELo^Snl" m the latter equations, the 
eigenvalues £ of the model (8) can be obtained as their linear terms in r\: 



2s+n 



£=£ (X)+J2 



2s+n 

- Xx o J2 t— "T" ~ A ( A ~ - Aa a + sa; 



(21) 



where £ (X) = (A - A)(Ay+ + sy ) - Xx^ + sx+; with x° = fig " ai *o5 ~ #g and Vn = ^g + ai? S a = {+> ~>- 
The quantities A& are fixed by 

^--(A fc -A)yo =«o E a^T ^ + ^ ' (22) 

where fc = 1, . . . , 2s + n; the quantum number n labels the excitations as discussed before. For generic 2s + n the 
equations above can be solved numerically. Alternatively the quantities Afc can be obtained as roots of the polynomial 
= IIsl~i"( A — A «) satisfying 

AP"(A) + — [y-A 2 + (y-A + x+)X - sx^] P'(X) - {( - (2s + n)X] P(X) = (23) 
x 

where £ is fixed by imposing that A = is a simple root of AP(A): ( = sx^ [17]. 

Conclusions. By the Quantum Inverse scattering method we have constructed integrable T-C models with twisted 
boundary conditions. The twist matrices are generic in the sense that they contain both diagonal and non-diagonal 
entries. They are responsible for the presence of rotating and counter-rotating terms in the Hamiltonian. The spectrum 
is computed through the Baxter method. As far as we know this method is applied to spin-boson systems for the first 
time; the subtleties related to the bosonic limit, recovered for infinite spin length are dug out. Integrability and exact 
solution can be obtained provided that a further term oc (a + a j( )S z is considered. Interestingly enough we found a 
global rotation of the spin/bosonic degrees of freedom such that the rotating terms (alternatively, the counter- rotating 
terms) are compensated out [29]. We conjecture that "true" counter-rotating terms in the Tavis-Cummings model 
could be inserted considering dynamical boundaries: Kx = Kx(X); alternatively one should consider XYZ symmetry 
of the scattering matrix. These terms serve to a reliable description of certain systems in quantum optics [26,27] or 
to model the spin-orbit interaction in heterostructures where the simultaneous Rashba and Dresselhaus terms are 
important [30]. Our paper could pave the way to costruct integrable Hamiltonians for such physical situations. As 
immediate application, we notice that the Hamiltonian (8) describes the quantum circuit of Fig.(l). Two coupled dc- 
Supcrconducting Quantum Interference Devices (SQUIDs) are coupled inductively. The primary device p is intended 
built with large Josephson junctions to be described by a classical SQUID Hamiltonian [31] (whose degrees of freedom 
are, then, bosonic) flowed by the current I p (this circuit plays the role of the LC resonant circuit of Ref. [11]); the 
secondary SQUID s, with small junctions, is accommodated inside the primary and pierced by the magnetic flux: 
$ = 4>ext + Lplp (L p is the inductance of the circuit). Thus the secondary is a quantum SQUID controlled by the 
classical one. 



4 




FIG. 1. The quantum circuit described by the Hamiltonian (24). 
flux-qubit s by the inductive coupling caused by $ = <f) ex t + L P I P . 



The primary device p is a resonant circuit controlling the 



The effective Josephson coupling of the quantum SQUID depends on the flux Ej($>) ~ Ej(<fi ext ) + L p I p . This 
kind of setups are intensively studied as controllable flux-qubits [32,33] to data-bus transferring in many protocols of 
quantum computation [11]. The circuit Hamiltonian is 

Hdrcuu = u p a^a - 2E 6 j(<p ext )S x - 2L p (a + d)S x - iM(a - a^)S y + 2V c (a + a))S z (24) 

where tu p is the "frequency" of the primary SQUID (or the natural frequency of the resonant circuit [11]), Vc is due 
to the capacitive coupling between the SQUID's; EJ, E c ,a = {p,s} are related respectively to the Josephson and 
the charging energies of the junctions and M is the mutual inductance; the gate voltage V g is tuned to the charge 
degeneracy point [31]. For generic circuit-parameters the dynamics of the qubit is intricated by the presence of the 
counter-rotating terms, making the device not reliable in the communication protocols. Our calculation suggests how 
the circuit parameters can be tuned to reproduce our model (8) ; for it the dynamics is not altered by the presence of 
the counter-rotating terms. Using this trick the qubit dynamics can be effectively "protected" at any frequency u> p . 
The relation between the circuit-parameters and coefficients in the Hamiltonian (8) is: {ui p , Ej, 2L p , M, Vc} — > {X(U+ 

V),2Y, U+V, U-V, -2VUV} (Z = is set for simplicity), implying that Vc should be tuned to Vc = \J M 2 - AL P /2 
to make un-effective the counter-rotating terms. 
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